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APPROVAL 


The contents of this documentation are in accordance with 
good engineering design and computer programming practices. 
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This documentation is correct and complete and the program is 
coordinated with related work in other branches. 
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“LIMITATION OF LIABILITY 


The U. S. Government will incur no liability resulting from use. 
of this computer program by organizations or persons other than the 
U. S. Coast Guard. 
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ales Vite: Jecbrveaher Propulsion System Simulation (Diescl-Flectric). | 
2. Brich Description: This program solves numcrically the set of simul- 

taneous non-lincur differential equations which describe the dynamic 4 

| 


behavior of the complete propulsion system (including vessel dynamics) 


of a diesel-electric icebrcaker during rnancuvering. ‘The program 


also contnins a means of subjecting the propeller to ice impact load- 


ings of the user's making. At the option of the user, input data can 


readily be prinied for decumemtation purposes and output data can be 


plotied, 


3a. Author: Jack W. Lewis, Licutenaiut, U.S. Coast Guar] 


b. Date: May 1968 


ec. Documentation: Prepared by Computer Usage Development Corporation, 


Washington Scientific Office. 


4, Language: FORTRAN 1V (IBM 1130 version) 


3. Machine: IBM 1130 with d’sk (plotter desirable but not necessary). 


6. Security Classification: UNCLASSIFIED 


hy astimated Running Time: Will vary directly with length of maneuver 


and time step employed. Using time steps of .08 seconds will require 


approximately 7 minvtes computation time for 15 seconds of ship time. 


After initial propulsion plant transients have disuppeared, generally 


within 10 to 20 seconds after the start of a maneuver, the time step 


can be increased to .15 to . 20 seconds thus reducing computation 


time. 


PART Il - PURPOSE AND METHOD 


Purpose: The purpose of this program is to simulate transient con- 
ditions within a diesel] electric propulsion system for an icebreaker. 
‘The program will simulate various conditions to which an icebreaker 
may be subject, including crash reversals and accelerations under 
free route conditions, ice ramming, and subjection of propellers to 
ice encounters. Use of the program will permit designers to real- 
istically study the effect of the transient conditions on design para- 
meters and to optimize design seleciion of equipment for new ice- 


breakers as rcquired by the U. S. Coast Guard. 
Method and Theory: 


a. Method: The construction of mathematical models for complete 
ship propulsion systums involves writing and solving a large 
number «f coupled non-linear differential equations. This com- 
puter model for sulving these equations describes as accuarately 
as possible the complete dynarnics of all components of the ice- 
breaker‘'s propulsion system. Numerical integration techniques 
are used to simultaneously solve this complete set of equations 
and to present the resulting data as a time history uf key variables 


both on a printer and graphic plotter. 


b. ‘Theory: The theoretical approach used in developing the model 
was to derive the equations and program them for computer 
solution. The computer solution for the USCGC GLACIER was 
compared to sea trial data for this ship. Complete theory and 
results of initial test runs is presented in references 1, 2 and 


3 below. 


References: 

ide ene Rat ak | 

ICEBREAKER PROPULSION SYSTEM SIMULATICN, U.S. ) 

Coast Guard Publication, 1967, by Lt. J. Lewis, U.S. Coast ' 

Guard, I. Scoville, and LCDR Ii, J. lecourt, U.S. Coast 

Guard. 

1) Phase I - Simulation of USCGC GLACIER. 

2) Phase IT - Dynamics of AC-DC Icebreaker Propulsion | 
System, | 


3) Phase Il] - Propulsion System Behavior During Propeller - 


Iee Interactions. 
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General Restrictions: See other restrictions. 


Tapes Required: None 


Nonstoudard Hardware: A Calcomp plotter is desirable but not 


necessary to run a simulation. Ifa plotter is not available, the 
GRAPH mainline routine and its four subroutines -GRID, DATA, 
VAPLO, and PRINT - must not be compiled nor any atlempis be 
made to execute the GRAPH routine (deta swilch 2 must there- 


fore reinain O!}). 


Maximum Array Sizes: The only size restriction of any concern 


¥ the Viinit of 600 records that can be stered on the disk during 


corun. A disk wrile interval tame step must therefore be such 


as to not cause more than GOO records to be written on the disk. 
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PART IV - NONSTANDARD OPERATING INSTRUCTIONS 
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1, Special Operating Instructions: 


a. 


General: This program has been written in a manner 
which places the engineer in the computation loop. This 
type of operation was choosen in lieu of batch processing |] 
because of 1) the experimental nature of the problem fre- 
quently requires decision making during runs and 2) the IBM 
1130 computer is quite adaptable to this type of operation. 
Considerable use is made of the keyboard console for the 
desicion making process. 

Documentation Procedures: The sole purpose of the mainline 
prosrei GNIS is documentation of each run, Although it is 
somewhat expensive in machine time to swap mainline programs 
in and out of core storage, the advantages of being able to docu- 


ment each run overides. 


Documentation is handled both in a decision making process 

at the keyboard console and by GENIS itself, In other words, 

certain input data and headers will automatically be printed 

as output wheras other input may or may not be printed depend- 

ing on the console selector switches turned on, On the very first 

run switch 9 should be turned on in order that the serial number 

disk file can be initialized, This switch should then be left in 

the off position until it is desired to reinitialize the serial 

number file. Printed instructions will appear on the keyboard - 


console to guide the user through each run, 


¥ 


Problem Solving: The primary problem involved in the 


solution of differential equations by numerical methods is 

the selection of the time step, If the time step is too large 
the numevical integration method ray lose the solution and 
"blow up.'' If the step is too small solution time can be 
prohibitive, Asa general rule it is best to make several runs 
at first with different time steps to compare accuracy and 
solution tine, Duxing the first parts of any mancuver (usually 
the first 10 to 20 scconds) experience with this program in- 


dicates time steps around ,06 to .08 seconds to be about right. 


Jt is possible to tnterupt the solution of a prohlem at any time 
by turning on console sclector switch zero. In addition to 
stopping the process of integration, turning this switch on 
calls in subroutine FINIS which sNows the user to select a 
number of optious concerning what he wants to do next, For 
example if it is desired to change the time step, switches 3 
and 5 would be turned on (make sure switch zero is off) and 


the push button marked "START" depressed. A message 


aK ) TIME STEP" will be typed at which time the vser 

types in the new value of the time step. The key marked 

"EO" (end of field) is depressed and the program continues i 
the problem solution using the updated value of the time step. i 


Experience has shown that the tirne constant associated with 
the armature circuit governs the size of the time step. Gen- 
erally speaking, after the armature circuit current stops 
fluctuating greatly (usually after 10 to 20 seconds of a man- 
euver) the tine constant of this circuit can be increased 
without introduction of errors, Increasing the time con- 


stant allows larger time steps and hence decreases solution 


o / 
| / 
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tirnc, Increasing this time constant can be readily 
accomplished by increasing the value of the armature 
circuit inductance (XA). <A value of .008 henries 
allowed a solution to be accomplished with time steps of 

. 20 seconds in one case with no noticable difference in 
the solution optained by not changing the inductance (, 0002 


hencries) and heeping the time step at .08 seconds, 


It is a simple programming problem to change subroutine 
FINIS to include additional options desired by the user . 


Switches 8 through 15 are still available for this purpose. 


Subroutine CNTRI. was written with the idea that the user 
could adda library of mancuvers and data changes through 
the use of the variable ICNTR. For example a "Flank 
Powcr Ahcad to Flank Power Astern" maneuver could be 
called using ICNTR equal to 2 and setting predetermined 
initial conditions is.to the program in addition to the ordered 


bridge controller position. 


It is often easier to obtain the necessary initial conditions 
for a particular maneuver by using the simulation. For 
example, suppose it is desired to conduct a "Flank Power 
Ahead to Flank Power Astern" maneuver. The initial 
conditions at "Flank Power Ahead" are needed. These 
could be calculated roughly by hand, fed into the simula- 
tion and then by keeping the simulation at "Flank Power 
Ahead" allow the simulator to find the correct initial con- 
ditions. Once found, a new set of initial condition cards 


can be created by using switch 1 option of FINIS, 


Program Alterations: Altcring the program has already 

been briefly mentioned above and will be further expanded 
upon here, An attempt was made to place components that 
may possibly be changed by a user into separate subroutines. 
This should facilitate changing components from one design 

to another, Jor example to change diesel engines, subroutine 
QDEV would have to be changed, If idle speeds or max speeds 
are changed it would also be necessary to make changes in’ 
subroutine CON2 which is the bridge controller command 


generator subroutine. 


It is hoped that the program contains sufficient comments 
so that a uscr with the aid of this documentation can make 
major or minor changes to the program to solve his specific 


problem 


Plotter Output: Mainline program GRAPH is called into 
operation by turning on selector switch 2 after a program 
interrupt. Plotting and labeling the solution data is the pur- 
pose of this program, Before the data is plotted the user 
must supply the program with three essentia] pieces of 
information, i.e., 1) which variables the user wishes to 


be plotted 2) whether or not a grid and labeling is to be con- 


structed and 3) any titling information. Data on variables to 


be plotted can be entered in any order; however, leave no 
spaces and insure the data is right justified. Drawing the grid 
requires time; therefore, during parametric investigation 

it is often best to first draw one grid and then use this for the 
remaining problems, Titling information may be omitted by 
depressing the key "EOF" three times. After the data is plotted 


mainline program GENIS is reloaded automatically. 
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; 
2. Mrvor Cosvectiows: None } 
3, Restart Instructions: None : 
4, Program Changes: Hach of the subroutines are written in 
toe eee ee 4 
= modular form and in FORTRAN, ‘This is done to allow the 
; user options of changing the routines to simulate special 
‘ situations, As an example of the type of changes that might 


be made let us suppose the armature circuit looked as follows: 


ES 
ea 
(s) © 
[ it ] 
4 
ik 
This circuit would only allow unidirectional current flow. Follow- 


ing is a modified subroutine CALVR which would simulate this 


situation - i,e. if EC >) EG then curvature circuit = 0.0. 
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PART V- DA'TA PREPARATION 


Input Data: The following indicates in detail how each data 
card is preparcd and from where the data is gathered or cal- 


ulated. 


Mach page or section of the report discusses individual data 
cards, preccded by an example data card prepared for the 


simulation discussed in references 1, 2 and 3. 
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PROPELLER DATA 
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Data-Card J es 
“ ee 3 3 eae ie: 
FORMA‘ 15, bs, GIF10, 4 
ELEMEN'IS é 


ZB = Number of Propeller Blades (3 for USCGC GLACIER, 
4 for M-12 hull) 


DAR= Pr ao developed area ratio (non dimensiona)) 


ee Cs ie. Oa 
1,067 - 0.229 x P/D ; 


a= 


DAR = = 


where DAR = Developed area ratio 
EAR = Expandcd area ratio 


. P/D = Propeller Pitch=Diameter ratio, 


One or the other (DAR or EAR) will be available from propeller 


characteristics. 


PDR = Propeller pitch/diameter ratio (P/D or H/D) - Non dimensional 


available from propeller characteristics. 


£ 


D= Diameter of-propellerfrom tip to tip (Feet). ae 


EO* = Thickness of propeller blade at tip (Feet). , 


Ap: "a8 Fae cas: baie x .- - 
eR he: ‘3 - oe -* 
of. : perigee * 


* See figure following the next page for dimensional locations on propeller. 
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KEK = Thickness of propeller blade at reot (ect) 


DCT = Dian der of hub (eect), 
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Hen BB(1) J(141) B(1+1) BB(1+1) 
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CAD = - Added mass coefficient (non-dimensional) 
Use 8%. i 

METTA = _ Ship displacement (long tons) 
XKK - Number of propellers. — 12 -6 
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VK= ' Speed of ship (knots) 
EHP = , Effective horsepower (horsepower) 
WADED Wake deduction (1-W,,) (non dimensional) 
THDED = Thrust deduction (1-T) (non-dimensional) 


Each of the elements are read from the ship curves such as those 

following, VK spacing must be as given as it is assumed by the 

program, not read, Values for EHP, 1-w, and 1-t not given’ on the. curve can 
be found from straight line measurements back to origin for 1- t and: lew 
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and a square law for EHP. 


Note: ‘Written description continues'on page 53 
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ELEMENTS 
Each of these elements is fouhd from data on the equipment such 


as from the certification daia shown for the GLACIER. Note that 


values are on a per shaft basis. 


RM = Resistance of motor armature circuit including 


the compensating windings. (ohms) 
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Ze number of armature force conductors 
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data for this is read from the technical manual or 


design sheets for the equipment. 
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ELEMENTS 
Each of these elements is found from data on the equipment such as 
from the certification data shown for the GLACIER. Note that values 


are on a per shaft basis. 


RG = Resistance of generator armature circuit including 


series field and commutating windings. 
‘ TRNFG = _ Turns of generator shunt field winding/pole. 


ZGP = Z‘ = flux conversion factor (volts /weber-rpm)from 
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RFG = ,.......... resistance of generator shunt field (ohms). 
ae Ble ONE ET eo 3 $ 


at ea @ o Suqprensing “OG 2 o See “HU OT = 
prvryyrrrTyy yt yyy) ™ (Toate) a acre vem (4) 
“WC = oy 70 eg (2) OE ote + aging ey 6) 


endand oe mat Renstpel wn fs 
an aS and OUD ee an Sse ont oa rere pee om omn*CUs <= ereqtg wee (8) 
Pee ne ee, res = ew OBOE Primoeng pe (0) 2 on 


law sa DOU © nas ort “otey onte “our ons Sn est (0 
“5US © ae eR “ete cag “aren fn “en ee Ce) 
GAR Sig » 


eels pee ey ete 


“Treg Pam yang 70 én St ews Spee am memme Ww 
Bepreds see Se RS INT pee “emtven srveease CoQ eey 


1 G 
4 i . | {LR SD ae an ma rat ae 
} irqet ee i aH : BREESE re sr : 
4 Ver t . ee ae . m rT T 
; oaseernee am: SaaREEn ana — Were t Se Fem Sree PRs ee ng Pt. 8 I (*) Ot we geaaituy “yas yameme om YO Fone = smeotna Fe ened (0) " 
RES ST ee Rae Tce RR MT. SEES a RS Sa ss 
ie ee : sis tm Ree eee vl 3 
: PRES Maret (imi ieee are kot $ 
Siaitatees Tos ee eae, 3 
| cents , : - ' 
; : , . - - ae 
i 2 4 if eee, ct a eR . i ss 
fONt ti if A i: ; i . is. aL © . 
1 i ee: : 4 e eo. ‘=a er a 5 | i w anpe aa: gan nina os nesta ees 
' —- : oo 6 A ced Aen Peeainn cee aan, pumen Wa Te UBT canteens (0) 
H ape a Ky ame SLOT (7) 
. a ee Ren ee sm ems oUece (0) 5 
“aSk ome (2 = eomonee (0) swtee encaferea (5) : 
z Seerreend © exe 8 SoEG SE OTR E Tt > Eee (0) sama outa (*) 
: | 
| ¢ 
' 
| 


-57- 


se 7 
nes 61 E wLOSeS Foret wees veTVe 
L ict . J ws ie — 
©. SEER setive SOC 22 Yew ou = yew Treat ce ea, (2) Sot = wry (0) Weare r a Prien oy riya 
"SC > et ery t2) “8 bo eg (2) Lac eaud < 
on wer -  WWULL ORS embes 250 0 CURLS es 09 JT OW 


Ceves CER wets 42RTT RE TU oqaaey Peg oye 


2 > oe et coe FP ewe (: 


SE cose = Se ary Toren (: 


= emg a) Ue we A * ee SB OTN" (5) 


8 cnr"e = Ge on Toes (YQ 


ee 4%. 14% 1eo5, De Vast WV as 


KIMAX Cl XK1 CSHFR TRNSE ‘RX 
CUOLEROO CECCOLOO OOCOEACOOSCONLOLCSOFOOHOLOCOOONHOCOOOOHGORCENNOOHOOVOTS O ba 
eee Seer erent renner herr ran rr ener rr acura sick ee Ree eA) 
PUTED TET GVTETE TEED A TTT ANTE TTD PETTITT DTT) of ott 


QLQVAULEUULTLDDDADLIDLPADAL GADD Z2QGAVDLADGAADADAAVDADADAADLDLZAALIAIA2 DVI 


™ 
™ 
N 
aa 
nO 
nm 


23553333 $39533539 363939599 3932393933 33333335 32233333 3335395333323593 3333 | 


3 
i CA ee ee ee ee ee 44 


_ 
- 
> 
~ 
xy 


CECE CCPECSECC COE RE GGOT OCC OLS COGCOCED COGCCOCCEECOCCOLGHCOESSOCHOCESCCCOEGECOCOCS 
CTUT UT TTT TTT ATT 00070777 oD 
CEGELEL F ESLELSCES GEHGRES BCEREEEREE CECLEEEED LALORER BELBTEREREECEGE SERA EDA 


GLEGEOLELEUSSGEUESTESGSESUSSSESLSGTS SOS UESSU TOTES STU NSCC ESTOS OSS I OSE IOS 
2I4SSCVO CUNY a UES SE Be ZRG o niet Peay eS SHER GS os Ha SS, SEVEESD TE GES CPCS LS ID G2 87 ED DS Bed 
“£OC £084 
Data Card 73" . 
FORMAT 61°10. 4 | 
ELEMENTS 
FIMAX = Strong field (bollard conditions) value of motor 
field current (amps) 
C1 = offset of motor controller curve (volts)(see - 
following figure) 
XK1 = slope of motor controller curve (amps/volt) (see 
following figure) ° “ 
CSHEFR = shaft frictional coefficient (foot-1b-sec”) . 
, 
ia) 4 
Gost © ~SHFR $58 ay 
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“TRNSF = generator series field turns/pole 4 
° 3 * rl 
“= resistance (ohms) of the vee iaile resistor in the 


generator field circuit. 
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On the GLACIER, the motor field is plugged in the reverse direction 


during a crash back approximately as follows: 


motor ; 
field ete ae eS TEN ae ER ee FMIRV + FMI, 
current ; . 


FMIRV = value in amperes per above sketch 


TAUA, TAUB = exponent values to fit the following equation to the above 
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The clements are derived from the generator and motor flux curves 


(see following figures) 


CFGA = constant times the generator field magnetic flux 
(volts / rpm) 
‘ CEMA = constant times the motor field ie uncles flux 
(volis/rmp) 
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“porMAT 0” aE 104 
ELEMENTS 
I : S. 
ALP1 = 4 buffer spring scale * 
ALP, = 1 Wy governor to engine constant* 
‘ ALP, = Wy needle value adjustrrent* . 
ee a ALPD = Wy .*. speed droop coefficient* 
XK6 = _ _ rated diesel load’ — Ss) s 4 
| . rated diesel stroke inch» seni 
The rated diesel load is read from curves such , oo 
as the attached. The rated diesel stroke is | 
: tL: | 
_ obtained from the manufacturer. 


* See reference 1 for discussion of diesel governor datz. 
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to friction 


QBRG 


see sketch bclow 


see sketch below 


a ee 


- 


— CWM=slope of line 
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Data Card 91. 
* “FORMAT TF 16,-4)" 
ELEMENTS 
XJPA = Polar moment of inertia of propeller in air 
2 
(ft-lb-sec ). 
If. given a value of zero on data card it 
will be calculated as per subroutine POLM I. an 
Otherwise, given value will be used. 
XJIM = Polar moment of inertia of one motor = 
) oe inom 2 
wk (ft - lb-sec ) : : a 
g- : uy 
where g = value of gravity i 
wk? = given inertia factor for motor 3 
Oe Nas 2 
# ' . XJSH = Propeller shaft inertia (ft-lb-sec ). Usually 
; . negligible, rae 
ee 2 ' 
XJID = Polar moment of inertia of one diesel (ft- lb-sec ). 
from manufacturers data. 
“KIG's 7 Polar moment of inertia of one generator = 
ear ; y) 
— (ft-lb-sec’) © . : 
g . ‘ 
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ELEMENTS 
RE = Rox = depth of cut factor for propeller 
D - encountering ice. Values are non 
: - dimentional for each successive pro- 
peller blade. See following figure. . 
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FGI = aA initial generator field current (amps) 
XNPM = initial diesel aoe (rps) 
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ELEMENTS 
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DT = time step (seconds) for integration 

'SIMSC = time increment (seconds) for writing solution 
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TICE = time (seconds) when start hitting ice 
oe 4 


TEND = time (seccnds) when program ends. 
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APPENDIX 


SYMBOL TABLE 


MATH 
FORTRAN SYMBOL 
LABEL _(if any) _* _DESCRIPTION 
RO p “Mass density of salt water equal to 1. 99 
1 2351 1 Constant equal to 3, 14159 ia 
IZB xf Number of propeller bladeg | | 
DAR DAR Fj Developed area ratio of propeller 
PDR P/)D ah Pitch-diameter ratio of propeller | | 
D 1») af Propcller diameter | | 
XJPA pair J Propeller inertia in air | | 
XJM ee af: Motor inertia | 
XJSH ah al Shaft inertia | 
XID y Diesel inertia 
XJG oS Generator inertia 
EFG Pee rdered generator field voltage 
XJIMP dip Inertia of propeller-motor system 
XJ DG Jomg - Combined diescl-generator inertia 
FIMAX ipmax J Maximum value of motor field current : 
RA R, Armature circuit resistance 
RG Re Jf Generator armature resistance | | 
ei C, af | gees used in calculating ordered | 
XK1 K, Ni armature current in motor regulator | | 
QBRG wv Y-intercept of developed diesel torque 
vs. angular velocity of prime mover 
curve 
CSHIFR Cy if Shaft friction constant 
CWM Pi Slope of developed diesel torque vs. 
angular velocity of prime mover 
curve 
PCMS af Entrained water inertia 


**a check(/) implies the variable is initially read in by subroutine REDAT. 
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FORTRAN 
LABEL 


XLA 
TAU1 


TAU2 


NOVAR 


NVAR 


JJ 
ICNTR 


NPLOT 


bce 4 


MATH 
SYMBOL 
(if any) 


L 
pi 


To. 


maw 


=e SS 
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ay 
A EA) i te vl 
= sna tev cet ota cond NaS ea cng eases, Aiea nimm eed Mt 
i 


DESCRIPTION 


Armature circuit ind:ictance 
Time for bridge coutruller to reach 


position 8, from initial position 6 ie 


Time for bridge controller to go from 


position 6, to the final (ordered) 


1 
position 8. 

Developed diesel torque constant 

Buffer spring scale 

Governor to engine constant 

Needle valve adjustment 

Speed droop coefficient 


Used to contain the value 2240. (1. +. ) 


4/32. 18 used in the computation : 
of the ship acceleration. | 
(Graph routine) indicates number of 
variables to be plotted 
(Graph routine) 10-word array con- 
taining the code numbers of up to 


10 variables to be plotted 


_ (Graph routine) index to NVAR array 


Indicates particular command given; 
equal to one (1) for ahead flank, 2 
for back flank, 3 for ahead full, 
and 4 for ahead full in ice. 

(Graph routine) indicates if a grid is 

: to be plotted; equal to zero (0) if 
a grid not desired, equal to one (1) 


if a grid is requested. 


a check (/) implies the variable is initially read in by subroutine REDAT. 


FORTRAN 
LABEL 


NDATA 
NON 


sYl 


SY2 


SY3 


SY4 


SY5 


SY6 


SY7 


SY8 


sy9 


MATH 
SYMBOL 


__(if any) 


DESCRIPTION 


(Not used) 

Indicates the number of time intervals 
between listings of the ovtput on 
the printer. 

(Graph routine) index to NVAR array 

(Graph routine) plotter scale factor for 
time axis | 

(Graph routine) pio:ter scale factor for 
variable currenily being plotted. 

(Graph routine) plotter scale factor for 
armature circuit current 

(Graph routine) plotter scale factor 
for bus volts 

(Graph routine) plotter scale factor for 
diesel speed 

(Graph routine) plotter scale factor for 
motor field current 

(Graph routine) plotter scale factor for 
propeller speed 

(Graph routine) plotter scale factor for 
power 

(Graph routine) plotter scale factor for 
vessel specd 

(Graph routine) plotter scale factor for 


propeller torque 


(Graph routine) plotter scale factor for 


propeller thrust. 


** a chock (,/) implies the variable is initially read in by subroutine REDAT. 
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MATH 
FORTRAN SYMBOL 
LABEL (if any) akel DESCRIPTION 


SY10 a . (Graph routine ) plotter scale factor for 
ship's distance 

SY11 m4 (Graph routine) plotter scale factor for 
ice torque y 

SY12 af (Graph routine - not used) , 

DT” > dt J Time interval | 

CT Save area for current time 

i. t af Current time 

TEND J Time ending value 

TIMSC J Time increment for writing values on 

the disk 

Z Z J Servo stroke position 

y¥ dz/dt Jf Servo stroke velocity 

XNPM pm f Angular velocity of prime mover 

FGI Jeg J Generator field current 

FIP i', J Motor field adjusting current 

XNP “as J Propeller angular velocity 

VS Ne A Ship speed 

XS vf Ship distance 

SIE Angular position of propeller 

Al is J Armature circuit current 

‘RV v (Not used) 

TEMP J (Not used) . 

CZ These 12 variables are related to the 

CY 12 variables immediately above 

CXNPM | "(relationship indicated by matching | 

CFGI labels after dropping the initial C) | 


** a check (,/) implies the variable is initially read in by subroutine REDAT. 
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| MATH 
FORTRAN SYMBOL 
LABEL (if any) 4% 
CFIP 
CXNP 
Cvs 
Ccxs 
| CSIE 
CAI 
CRV 
CTEMP 
: ZDO dz/dt 
| YDO ann fat” 
XNPMD ane "| 
| FGIDO di, g/t 
| FIPDO di! gi at 
f 3 . 
| 3 XNPDO dn / dt 
| | VSDO dv /dt 
| XSDO Vv, 
| | SIEDO n, 
[ AIDO di, /dt os 
RVDO 
TEMDO 
XNE n 
e 
XNEDO dn, /at 
VA V 
a 
THETO ox 
THET1 A of 
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DESCRIPTION 


(continued from previous page) 
and are used to save the values 
of the above variables at each 
step when the next set of values 
are being calculated. 
(CRV and CTEMP are not used) 


Servo stroke velocity 
Servo stroke acceleration 


Tine rate nf prime mover speed 


‘Time rate of generator field current 


Time rate of motor field adjusting cur- 
rent 

Propeller acceleration 

Ship's acceleration 

Ship's velocity 

Propeller angular velocity 

Time rate of armature circuit current 

(Not used) 

(Not used) 

Diesel error speed 

time rate of diesel error speed 

Advance velocity of propeller 

Ordered (final) bridge lever position 

Bridge lever position reached in? 


seconds 


** a check (,/) implies the variable was initially read in by subroutine REDAT. 


MATH 
FORTRAN SYMBOL 
_LABEL _(if any) = 
THETA @ 
THETD de/dt 
THEIC i, f 
Qne Cae 
G Q 
8 g 
EG e 
g 
EB e. 
EC e 
c 
QM Qua 
QFRM Seem 
AIHAT _ 
P 
Q on 
QIP 
QI 
! TP T 
p 
XPROI T 
| R R 
VAR 
XKW 
| CVK 
CFG c 
. fg 
| SPGEN if 
| CFGA ‘f 
TRNFG Neg J 
ZGP cae of 


* a check (,/) implies the variable was initially read in by subroutine REDAT. 


Fe eT ee ee ee ee ee ree) 
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DESCRIPTION 


Actual bridge level position 


Time rate of bridge lever position 


Initial value of bridge lever position 


Developed diesel torque 

Generator torque 

Generator vcltage 

Bus voltage 

Back emf 

Motor torque 

Shaft friction torque 

Ordered armature circuit current 

Propeller torque 

Ideal propeller torque 

Ice torque 

Total propeller thrust 

Thrust per prcpeller 

Ship's resistance 

(Graph routine) contains variable 
retrieved from the disk 

Ship's power in kilowatts 

Ship's speed in knots 


Generator flux 


Spacing of generator, independent axis 


for CFCA array 


13-word array containing generator 


4 
i 


ig 


flux for spacing indicated by SPGEN 


Number of turns in generator 


Generator flux cunversion factor 


re 
LL a, STP 


FORTRAN 
LABEL 


RFG 
ALNG 
TAUIFG 
TRNSIF 
RX 

XM 


XILAS! 
AMP 
AG IN 
CIleM 
SPMOT 


ClMA 


TRNEM 
ZMP 
RFM 
XLIM 
TAUFM 
FMIRV 


TAUA 
TAUB 


B 


BR 


MATT 
SYMBOL 
(if_any)_ 


Het 


nN 


om 


nN 


J 


-108- 


Generator field resistance 

Generator field inductance 

Generator time constant 

Number of turns in generator field 

Variable resistance in generator field 

Mutual inductance between generator 
separate and in serics fields 

Inductance of generator scries field 

Ampere turns of generator 

Number of generators 

Motor field flux 

Spacing of motor, independent axis for 
CIMA array 

13-word array containiug motor flux 
for spacing indicated by SPMOT 

Number of turns in motor 

Motor flux conversion constant 

Motor field resistance 

Motor ficld inductance 

Motor field time constant 

Field current increment for reverse 
direction 

Time constants controlling how fast 
IK MIRV gets added to FMI 

64-word array containing propeller 
thrust coefficients using J 

64-word array containing propeller 


torque coefficients using J 


MATH 
FORTRAN SYMBOL 
_LABEL (if any) hue 
E 
oe af 
El Cc 
Q af 
EHP 


WDED (1-w) 
WADED f 


TDED (1-t) 
THD.) A 


XkK K Jf 


TANT 
¢ Jan 


JEL 
VERIO 
a ©, 
MELLO 
SPARE 
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DESCRIPTION 


32-word array containing propeller 
thrust coefficients using 1/J 

32-word array containing propeller 
torque coefficients using 1/J 

20-word array containing effective 
horsepower for hull 

Wake deduction for hull 

20-word array of wake deduction 
factors for hull 

Thrust deduction for hull 

20-word array of thrust deduction» 
factors for hull 


Number of propeller shafts 


.Disk file record index 


Controls initial computation of ice torque 
Indicates position of ship in ice flow 
(Not used) 
(Not used) 
5-word array, the first 3 words of which 
are used as follows: 
SPARE(i) equivalenced to FMI (inn? 
-motor field current 
SPARE(2) equivalenced to ANGR(n,,) 
-ordered diesel speed 
SPARE(3) equivalenced to XNGRD 


(dn, fat) -ordered diesel speed rate 


** a check (V) implies the variable was initially read in by subroutine REDAT. 


on 


MA‘ 
FORTRAN SYMROI, 


_LABWI, — _(if any) = DESCRIPTION | 
ISPAR 5-word array, the first 2 words of 


which are used as follows: 
ISPAR(1)-controls calling of GENIS 
mainline from lXENC mainline 
ISPAR(2) -controls motor field 


plugging in CALVR 


CAD Cad a "Added mass" coeffieicnt of hull 

DUTTA A J Vessel displacement 

kM xi Motor armature resistance 

RBA ff 29-word array containing random 
deptlis of cut propeller blade makes . 
in the ice as a percentage of the | 
diameter | 

RBMX Rex Maximum depth of cut the propeller blade 
makes in the ice 

BLOC1 J Length o first block of ice | 

DBBI12 Bi Distance between ice blocks 1 and 2 

BILOC2 MA Length of ice block 2 

DBB21 af Distance between ice blocks 2 and 1 in 
repeating sequence 

FO ey af Tip thickness of propeller blade 

hk or J Root thickness of propeler blade 

BT Det J Diameter of propeller hub 

DCR P J Crushing pressure of ice 

DSH Do af Shear pressure of ice 

CICH af Corstant used to increase hull resist- 


ance for ice conditions. 


* a check (/) implies the variable was initially read in by subroutine REDAT. 


roo) 


MATH 
FORTRAN SYMBOI, 


_LABEL (if any) we DESCRIPTION 
TICE J Relative time from beginning of run vntil 


ice ig encountered 


CK1 J (Not used) 
CK2 J (Not used) 
CK3 J (Not used) 


** a check (/) implies the variable initially read in by subroutine REDAT. 


et 


SYSTEM FLOW CHARTS 


The following flow charts outline the flow of control in 
the three mainline programs (GIENIS, EXEC, and GRAPII) 


during a simulation run. 


a ieee ee tt inmate Dina 


Sait saan anaes 


GRAPH MAINLINE 


GRAPH 


Read In Variables 
To Be Plotted 


Subroutine GRID 


| Plot Grid with 
_,,! titles and scale 
information 


“Subrouti ne VAPILO GENIS ) Subroutine 


“bi Retrieve scale REDAT 

factor for vari- Read All 

able to be plotted Data 
paca . Cards 


GENIS MAINLINE 


Subroutine DATA 


Yetrieve valuc of 
variable from 


Subroutine 
' POLMI 


—— 


Inertias not 
Read In 


Plot Variable | 


a 


Subroutine PRINT 


and Requeste 
Input Data 


EXIT VIA LINK 
TO EXEC 


EXIT VIA 
LINW TOGENIS 


EXEC MAINLINE 


Subroutine 
FINIS 

Process 

Requests 


LXTT VIA 
LINK TO 
GENIS 


EXEC 


Write on dis 


|? 


a ae nee 


“ae Print OUT 
—— current 
valucs 
of variables 


Subroutine CNTRL 


Determine type 
of mancuver 


Subroutine 
CALVR* 
Calculaic 
initial values 
| of variables 


Save initial values} 


Subroutine OUT 
Print OUT 
initial 
values 


Subroutine ane Subroutine CALVR* 

pdate time step Calculate values 

and caluclate new of variables for 
updated time step 


values of variablesf* 


(called 4 times for 
each time step) 


Subroutine OUT 


So ee 


* Detailed in a separate flowchart 


SUBROUTINE CALVR ; 


( calvR ) ; 
ig 


Subroutine FILUXM 


Subroutine CON 1 
Calculate actual 
and time rate of 
bridge lever 
position 


Calculate motor 
field FLUX, time 
Hconstant, and 
inductance 


in ee 


Calculate 
motor field 
current 


ee nega scene 


Subroutine CON2 


Calculate ordered 
generator field 
voltage, dicsel 
speed, and dicse] 
speed rate 


Subroutine HUI.L 
Calculate ship's 
resistance, wake 
and thrust 
Ideductions 


Se see ae 


Calculate 
dieselerror 
speedand 


Calculate 
advance ye- 


locity of pro- 
peller 


Nea 


Subroutine QDEV tine rate Subroutine PROP 


Calculate develope c 
diesel torque 


Calculate ideal 

propeller torque 

jand thrust per 
propeller 


Adjust devel- 
opcddicsel 
torque for 
friction torqu¢ 


Calculate 
total propel- 
ler thrust 


Subroutine QICE 


—: 


Calculate ice 
‘torque of the 
propeller 


Subroutine FLUXG 
ae 
Calculate generator, 
FLUX, inductance, 
and time constant 


Subroutine |.RBMAX 
Calculate inax. 
depth of propeller 
cut in ice 


Subroutine STATE ee 

Calculate Derivatives: 

1. SERVO Stroke velocity 

2. SERVO siroke acceleration | 

3. Time rate of prime mover speed 

. Time rate of gencrator field current 
Time rate of motor field adjusting 

current 

}. Propeller acceleration 

. Ship's acceleration 

Ship's velocity 

Propeller angular velocity 

. Time rate of armature circuit current 


ee ES 


MAINLINE PROGRAMS AND SUBROUTINE 
FLOW CHARTS 


The following flowcharts outline in detail the processing associated 


with each mainline program and each subroutine. 


GENIS MAINLINE 


GENIS 


Print out 

peel con- 

eis We 
data 


Call DATSW 
4, KEY4) 


REDAT. 


Inurtias 


Print Title ] out 

Page awitch set- 
ting in- 

structions 


Call DA TSW 
(9, KEY9) 


PAUSE 


Call DATSW 
{11,7 


Reset 
serial No. 


Update Last 
Serial No. 


and title 
page 


DP iicecenieinciaeend te nail 5 7" 
oF. 


inkl ‘NaS ¥ eh eae 


GENIS MAINLINE (continued) 


Call DATSW 
(42, J) 


Call DATSW 
14, L) 


(13, KEY13) 


Call DATSW 


OFF 15,M 


Print out 
systems 
inertias 


Call DATSW 
10, KEY10) 


Print out 
Tce data 
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SUBROUTINE POLMI 


21.024 | D y5 
Pad 32, 18'1, 33 


*(1, OODAR 570.23 


ou 


pale 
1,33 


yt 024 (~2.)°(1, 37DAR r.., 30) 


Jap im sh’. Pair’? Pad 
J. =(XJD+XJG)XGEN 
pmg 


RETURN 


SUBROUTINE REDAT 


REDAT 


Read Prop _ 
Data 
(Cards 1-49) 


Read Hull 


Read Governor-, 
Diesel 
Data 
(Card 90) 


Read Systems 
Inertias Data 
(Card 91) 


Read Ice 
Conditions 
(Cards 92-97) 


ditions and con 
tro] data 
(Cards 


RA=RM+RG/XGEN 
XMASS#(14C, 3)2240+ 0/32. 18 


Seat ary 
RETURN ) 
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a ES Ee Ty ee 


KXEC MAINLINE - ee 
EXEC 


Initialize 
Control 
Variables 


CNIRL \ 
Compute 
Initial 
Values 


Call DATSW g| 


(O, KEYO)_ 


ON 


FINIS 
Clean up 
Routine 


OFF 


ISPAR(1) 


/_HEUN_\ 
Update 
Call LINK 
GENIS 


EXEC MAINLINE (continued) 


Convert Values t 
Proper Units 


Oe a a a ee 


SUBROU'TTLINIG CNTRI, 


CNR 


fomunanc 
Given 


‘ 


fame meee ee eter eet ee ee re eer res emer cnemraaranae seme = 9 


| a Se, 


eee Ny SS asl 


Ahead tank Back Plarile’ Ahead Full Ahead Full 
0,= 165 O05 -165° @, = 135° In Ice 
| 9,2 5° 


—CAIVR 
Compute 
Initial 


Valucs 
Save All 
Values 


Print Col. 
Headings 
For 
Output 
Data 


ey a 


Write Initial 
Values On 
Disk 


T 


( RETURN ) 
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SUBROUTINE HEUN 


HHUN 2 


set: 
DD=h/2 
Compute: 
sh: fm: 
Kk, h {(X y,) (4) 
Set: IJAX=1 


| t=CT+DD | 


CALVR 
Update 
Values 


Compute: 
ky=h-f(X +h/2, 
Save current values 
¥ tk, /2) of variables, 
Set: IJAX=2 Save tin CT. 


aoe 


&) 


Update 
Values 


Compute: RETURN | 
k,=h-f(X_+h/2, (4) 
3 n 
¥tk,/2 


IJAX=3 


et eee | 


}Set: 


_ i — i taille acre cari alia aaa tila i i i ali li ie ae a atari ee 
i aed 


SUBROUTINE HIEUN (continucd) 


| | t=CT+dt 


CALYR \ 


“Update 
Values 


Compute: 


Ya Yk tek 


2 
H2k,tk,)/6 
Set: IJAX=4_ 


[re |: 


SUBROUTINE CALVR es 2 — 
f__cont _\ 
COMPUT 

0 AND 
do/dt 


= CON 
COMPUTE 


7 


Ns dN. / dt 
BY er 
\ hen c 
aes 
COMPUT: © 
N =N. -N 
ec gr pm 


! > - J 
she aN, fd OM rat ot 


QDEV__ 


“COMPUTE 
Q 


de 


ee ememnen ane, 


CG TgaRGiewin WRN 
ye 7 ge QBRGFCWMEN XG) 
: i pyl 2- TRNSE i, (2: XGIEN) 


oe es 
AMPTR Nig 


SF 
COMPUTE 
GENERA TOR 


LUX 
VALUES 


COMPUTE: 
e =60'C. -N 
g fg pm 


FILUXM 
COMPUTE 
MOTOR FLUX 
VALUES 


a ama RR NT I TET aT RR ATO PTD 


SUBROUTINE CALVR (continued) 


COMPUTE: 
e,=60° Com’ Ny 
Qgr7. 04Cp," i 
Qm* =7,04C)_ +i, 
b*%g g AR, | SGEN)i, 


SET: 
TTN= t 

FMIIC=igy, 
IPLUG=2 


i, =FMICHI -T 


fm Ta" Tp 


os 


Stem Oe [18 


COMPUTE 
TDED AND 
WDED 


V_=WDED:V 
a s 


ds-- 


SUBROUTINE CALVR (continued) 


PROP 

COMPUTE 

QPI AND 
XPROI 

en 
-QPI 

2 Q 

'Y ="1DED-XPROI-K 
p w 


__ STATE 
COMPUTE 
DERIVATIVES 


SUBROUTINE CON] 


COMPUIS: 
@=0 +9, -@, ) 


IC ie 
14 
lao 1-10? 
dt 7 
COMPUTE: 
8-0, Ht-7,)(0%-0,) 
I 


RETURN 


| 4 
t ld 
14 
-129- 
4 
| oo 


SUBROUTINE CON2 


COMPUTE: 
=3- + 
34. UR, Rea? 


9 
tty 


max 
THET=|0| 


COMPUTE: 
ec, =e 


fg 'e ax 
(420. ppyyaepy¢ 29. * 139: 490, 9) 
ee ae 


gr. __420 0 
(72. x 60, jdt 


COMPUTE: 


e, =e 

fg fE ax 

go, JHET+(810. Bo 90. )) 
30. 


COMPUTE: 


9 
“tg THET fg 


RETURN 
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SUBROUTINE FLUXG 


.- FLUXG 


_/AMPTR| 
SPGEN | 


KGsIFIX(IF°G) 


FG 


Fe ah 


DCGDI=CFGA(13) 
-CFGA(12) 


DCGDI=CF'GA(KG+2)-CFGA(KG+1) 


oe CFGA(KG+1)+DCGDI: (FG-KG) 


——— ee: 


A a AMET G 
fg IAMPTRI “fg 
2 DEGuI _ 


XM=N,_- TRNSF-—2CG2E_ 


tg z," SPGEN 


Y.LASF=TRNSF*- —2CGD! 


Z_*SPGEN 
& 


RETURN 


— 


SUBROUTINE FLUXM 


FLUXM 


i 
| fa 


FM" opMOT 


KM=IFIX(}*M) 


Com CEMA(18) 


DCMDI=CFMA(13) 
-CPYMA(12) 1° 


DCMDI=CEFMA(KM+2)-CFMA(KM¢ 1) 
Con? CPMA(KM11)+DCMDI- (FM-KM) 


ZS BEML 
m 


fm F, Lom! Rem 


( RETURN 


SUBROUTINE QDEV 


=( 


14000: 0. 62 | ei 


1750 


=(———-60N xX GEN 
© max ( 90 ma 


Qe 


"Qe 
max ’ 


HQBRG+CWM-N,, )KGEN 


SUBROUTINE TLUEI, 


Vie | 
Vet) eo 
KV=V 


REHP=BHNPIRV41)4+(8 HP(KV+42)-T PCS V+ 1))(V-KV) 

WDED=W ADIUDUIRV+1)4+(WADEDUS V+2)-WADED(IKV+1))(V-KV) 
TDED= TMDED UK Vt1)+(THUDED(KV+2)- THDED(RV+1))(V-KV) 
REIJIP -550 

Ra---- cearees a2 (Sey see 

Vs 


R=1,4R 
WDED=1, 
TDED=, 9 


oe 
=o 
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N 
{ NYOLAY y NOON LUY 


a 


RS errr ce 


@0Od -2VA -‘ALYX=10OUdX 
d-2dOU-,VA -dOxX=1d 


os 
£0> 
a 
0<~ 5 
x 3 
fX-01T=ds 
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SUBROUTINE QICE 


} Le] 


/ REMAX 
Compute 


Rie 


COSla, 19) 1-(2' Ry /D) : 
SIN(a, )2)*(1-COS(a4, 7 9)) 
O72" ARCTAN(SIN(0, 9)/COSla 
ABB=2:1/Z, 


b/2” 


RETURN 
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SUBROUTINE QICE (continued) 


COS) COR(a, 15) COS OTE) HSIN G5): SIN(SIE) 
R =D/2(1-COS(a, ; 9)/COS(B)) 


b 
RT (D-R,)/2 


a -ARCTAN(V /(R_- a.m N,) 
a, =(P/D): D/(2n-R_ ) 
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Be oo ND [(22Z, (D-D_,)) 
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SUBROUTINE STATE STATE 


dre on. dz, 
—z = 600 gat 1 O0%) & 3 nto ite oar “QAO 


N 


Hed = oy + ° 
VOLT 1 ere (Reg Rigg 
VOLT2=c -e -R ‘i 
f G a a 
XLD=L +1 .5/XGEN: XLASF 


DENOM=1. 5 Lig LD-XLM 
di 1,5(XM: VOLTI+Lg,," VOL/T2) 
DENOM 


1.5 XM: VOLT24XLD-VOLT1 


dt XMASS 
dx, 

dt a : 
SU aN 2n 


RETURN 
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SUBROUTINE OUT 


RPMPM=60- N 
pm 
RPMP=60- N, 
CVK=V,/1. 889 
TPT=T 
p 


RT=R/2240 
POW:Q_2nN_/550 
mp 


PRINT OUT: 
t, RPMPM, leg? RPMP, 
CVK, Que’ Qe “> tm 


RT, POW. 


RETURN 
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SUBROUTINE FINIS 


FINIS .) 


Call DATSW 
(4, Ky 4) 


Call DA TSW 
_(&, KEYS) 


OFF (=2) 
ON(=1) 


ON(=1) 
OFF(=2) Read In 
—_—_ New dt 
tat ; 
17 


Ly “Call DATSW_ 
oath ieee ee 
| pausn (6, KE Y6) 


Call DATSW 
_UO, KEY) 


OFF (=2 


printing 


interval 
if (NON) 


PUNCH AND PRINT 
NEW INITIAL, VALUES: 


1.0.0 ON ny! tt dz/dt, z, 


Call DATSW 
(7, KBY7) 


; as 
‘fe’ Nam! *f? p 


Call DATSW 
(2, Ki Y2) 


Read In 
New L 
a 


Call DA TSW 
(3, KEY83) 
ISPAR(1)=KEY3 


“ ON(=1) 


| TEND=t-di | 
; awed “ae 


re enemne <aemeerel oenree 


aia = 


Call LINK RETURN 
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SUBROUTINE GRID 


( GRID 
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Pin 
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SUBROUTINE GRID (continued) 
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SUBROUTINIE GRID (continued) 


SIT POSI'TION 
FOR NEXT 
TITLE ON 
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RETURN PEN 
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UPDATE INDEX RETURN 
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SUBROUTINE VAP).O 
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SUBROUTINE PRINT PRINT 


Retrieve 
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RETURN 


MAINLINE PROGRAMS AND SUBROUTINE 


DOCUMENTA TION 


Program Name: Mainline Routine GENIS 


Program Description: The GENIS mainline routine, after 
calling subroutine REDAT to read in all the input data and 
subroutine POLMI to compute the remaining inertias, prints 
out the control data and initial conditions and depending on the 
setting of data switches 10 through 15 prints out the ice data, 
propeller data, hull data, electrical data, governor and bridge 
controller data, and inertias respectively. The EXEC main- 
line routine is then loaded and executed (by use of a CALL LINK 
(EXEC) statement) to begin the simulation. 


Program Logic: Disk data file 2 contains the serial number of 
the last simulation. If data switch 9 is on, the serial number will 
be reset to zero (switch 9 must be on for the first run), otherwise 


the serial number will be updated by one (1). 


Listing of Parameters: 


Values taken from COMMON - All values read in by subroutine 
REDAT plus XJPE, the entrained water inertia which was computed 
in POLMI. 

Values taken from the disk - SER, the previous serial number if 
data switch 9 is off. 

Results stored in COMMON- SER, the updated serial number 


corresponding to this run. 


Program Length: 


a. Source program - 245 statements (including comments). 
b. Object program - 1828 machine words. 
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6. Programs Called By This Program: Subroutines REDAT and 
POJ.M1; the EXEC mainline is loaded and executed through use 
of a CALL LINK (EXC) statement. 


he Programs Calling This Program: This mainline is loaded and 
executed cither by 1) the initial execution as specified on the IBM 1130 
Monitcr Control Record // XEQ GENIS; 2) the EXEC mainline upon 


the completion of a run; or 3) the GRAPH mainline upon the com- 


pletion of a graph. 


8. Local Variables: 


LIST - 40 word integer array used to store the title 
(80 characters or ieee) ‘ 

KEY4, KEYS, KEY10, KEY11, KEY12, KEY13, KEY14, and KEY15, | 
integer variables used to store the setting of the 
appropriate data switches; 1-equals on, 2 equals 
off. 


eo} 


Printed Output: The following lines are printed on the 1132 printer 
as a title page - . | 
UNITED STATES COAST GUARD ICEBREAKER PROPULSION | F 
SYSTEM SIMULATOR | d 
SIMULATION OF (80 character title) | 
| 


SERIAL NUMBER xxxxxxxxxx 
At the top of the next page is printed - SELECTED INPUT DATA 


followed by the following headings with their appropriate data- 


iNITIAL CONDITIONS . 
CONTROL DATA | | q 


PROP DATA 

HULL DATA 

ELECTRICAL DATA 
GOVERNOR AND DIESEL DA'TA 
SYSTEM INERTIAS 

ICE DA'TA 


Finally, at the top of a new page is printed OUTPUT DATA. 


The following instructions pertaining to the sclection of the 


output data is typed out - | 


TYPE IN CONCLUSION OF PHRASE 'SIMULATION OF' 
USING 80 CHARACTERS OR LESS. 


TURN ON APPROPRIATE SWITCHES 

ICE DATA -10 

PROP DATA-11 

HULL DATA-12 

ELECTRICAL. DATA-13 

GOVERNOR AND BRIDGE CONTROLLER DATA -14 
INERTIAS-15 


TURN ON SELECTOR SWITCHES AND PRESS START. 


1 Program Name: SUBROUTINE POLMI 


cs Program Description: Subroutine POLMI computes the inertia 
of the propeller-motor system Or and the combined inertia 

of the diesel engine and generator armature (J D: 

pmg | 4 

} 4 

3. Program Logic: The number of propeller blades (IZB) determines | 
which of two formulas will be used in computing the "added inertia" 
Pad” If a value for the propeller inertia | 
in air (Spain) was not entered in, its value is computed. The values of ‘ 


of the propeller in water (J 
J and J are then calculated. 
mp pmg | 


4, Listing of Parameters: All parameters are in COMMON. 
Values taken from COMMON - IZB,XJPA, XJM, XJSH, XJD, XJG, 
XGEN, and D. 

Results stored inCOMMON - XJPA, XJMP, XJDG, and XJPE. 


5. Program Length: . | q 
ie apse nee ne cae nae i 
a. Source program - 45 statements (including comments). 
Dd. Object program - 112 machine words. 


6. Programs Called By This Program: None 


re a rr re tS 


die Programs Calling This Program: GENIS mainline. 


8. Local Variables: None 


9, Printed Output: None 


cards that represent the propeller data, hull data, electrical 
data, governor-dicsel data, systems inertias, ice conditions, 


scale factors, initial conditions and control data. 


Program Logic: After reading in all the data cards, the program 
sets the values of p and 7 to 1.99 and 3, 14159 respectively. The 
armature circuit resistance is calculated from the resistance of the 
motor armature (Ray? the resistance of the generator armature 
(RQ), and the number of generators (KGEN). The "mass"term is 


also calculated. 


Listing of Parameters: All parameters are located in COMMON, 
Values taken from COMMON- none. 


"Results stored in COMMON - the variables initially read in by this 


program are so indicated under the COMMON data description section, 


Program Length: 


a. Source program - 127 statements (including comments). 


b. Object program - 498 machine words. 
Programs Called By This Program: None 
Programs Calling This Program: GENIS mainline. 


Local Variables: None 


_Printed Output: None 
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ee ee ee 


Program Name: Mainline Routine EXEC 


Program Description: The EXEC mainline routine coordinates 

the actual simulation. Subroutine CNTRL is called to compute 

the initial values afterwhich subroutine OUT is called to print 

them out, ‘hereafter subroutine HEUN is called each time a 

new time step is taken. Subroutine OUT is called each time the 
values are to be printed out. Upon conclusion of the simulation 
subroutine FINIS is called to select control for the remainder of the 
run. If FINIS returns(implies a graph was not desired) then the 
GENIS mainline may be loaded and executed if it is desired to 
start a new problem, otherwise HEUN is called to continue the 
present problem. The EXEC mainline writes the values on the disk 


for use by the GRAPH mainline. 


Program Logic: If data switch zero (0) is on after the initial conditions 
are computed subroutine FINIS will be called, otherwise it will 
not be. The variable NON controls the ship interval for printing with 


the variable MIKE being a counter of the steps taken. 


Listing of Parameters: All parameters are in COMMON . 
Values taken from COMMON - TIMSC, NON, CT, TEND, XNPM, 
XNP, CVS, QM, PI, AI, EB, FMI, CXS, QPI,XPROI, IELO,QI, and 
ISPAR(1). | 


Results stored in COMMON - JELO, JELO,KELO, and ISPAR(2). 


Program Length: 


a. Source program - 108 statements (including comments). 


b. Object program - 264 machine words. 


HEUN, and FINIS. The GENIS mainline may be loaded and exe~- 
cuted through use of a CALL LINK (GENIS) statement. 


Programs Calling This Program: Loaded and executed by the 


GENIS mainline through use of a CALL LINK(EXEC) statement. 


Local Variables: 


iI - initially set to one (1) but will take on a 
value of two(2) the first time the value of 
the current time (CT) exceeds 20 seconds; 
is used to reset the disk write time step (DTT) 


to 1.0 after CT exceeds 20 seconds. 


MIKE - initially set to one (1) but is updated by one 
each time a new time step is completed; is 
used to count time steps for comparison with 


prini out tisae step (NON), 


DTT - disk write time step; is mutially set to the 
value of TIMSC cut will take on the value of 


one (1) after CT exceeds 20.0 seconds. 


Tre = time value for next writing on disk; is 
updated by DDT each time the values are 


written on the disk. 


KEYO - contains setting of dita switch zero; 


equals one (1) if on, two (2) if off. 


RPMPM- angular velocity of prime mover in revolutions 


per minute. 


RPMP - propeller angular velocity in revolutions per 
minute. 

CVK - ship's speed in knots 

XKW - ship's horsepower 

JENE - set to value of ISPAR(1) for use in computed 


GO 'TO Statement. 


ie Printed Output: None 


-159- 


| 1; Program Name: SUBROUTINE CNTRL 


2. Program Description: Subroutine CNTRL sets the value of 


( the ordered bridge lever position (0) depending on the value 
of ICNTR which indicates the command that was given. Sub- 
( routine CALVR is then called to calculate the initial values, 


afterwhich the current values are saved. 


3. Program Logic: The value of 85 is set to 165.0 if ahead 
flank was given (ICNTR=1), 135.0 if either ahead full (ICNTR=3) 
or uhead full in ice (ICNTR=4) was given, and -165. 0 if back 
| flank (ICNTR=2) was given. The current values of the eleven 
ks 
pm’ We leg Mp V5 XS, and SIE) 


are placed into alternate locations, The initial values are then written 


state variables (t, ins dz/dt,z,n 


wey 


on the disk and a heading is printed out to label the columns for. 


succeeding print outs. 


[ 4, Listing of Parameters: All parameters are in COMMON. 


Values taken from COMMON - ICNTR, FIMAX, FIP, T, Al, Y, Z, 
[ XNPM, FIP, FGI, XNP, VS, XS, QM, PI, EB, PI, XPROI, QI and 
: IELO. 
| Results stored in COMMON -SIE, THETO, FMI, CT, CAI, cy, CZ, 
CXNPM, CFIP, CGI, CXNP, CUS, CXS, CSIE, and IELO. | 


[ 5. Program Length: 
( a. Source program - 84 statements (including comments). - 


b. Object program - 242 machine words. 


| 6. Programs Called By This Program: Subroutine CALVR 
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Programs Calling This Program: IEXEC mainline 


Local Variables: 


RPMPM.-. angular velocity of prime mover in revolutions 


per minute, 


RPM? - propeller angular velocity in revolutions per 
minute, 

CVK - ship's speed in knots 

XKW - ship's horsepower 


Printed Output: Prints out the column headings for the output 


data that appears in the printer as follows: 


T(SEC) RPMPM CFGI RPMP VK QDE(FT-LB) 
QG EB(V) FMI EC QM(?T-LB) AI(A) 
QP THRUST(T) RESIST HP /SH 


Program Name: SUBROUTINE HEUN 


Program Description: Subroutine WEUN is the numerical 
integration routine. The Runge-Kutta 4th order method is 
used. Very briefly, the method is as follows: we are seeking 
to solve the equation dy/dx = f(x,y). Given an (x4 Y,) we 
4s 
alculate: =y + + 2k. + +i 
calculate bare fe (k, 2 9 2k. k,)/6 
where: ky = af(x,, v0 
= -+- +k 

ky hf(x h/2,y,, k,/2) 

— 5 +h/s + 

ky hf(x h/2,y,, k,/2) 


kK, =hfi(x + + 
k, hf (x, h, ye ky) 


_ and h = independent axis increment (dx). 


a check is made on the values of Z and FIP. If Z is greater than 
1.000 Z is set to 1,000 if FIP is negative, FIP is set to zero. 
The local variable IJAX is used to return control to the point in 
the routine where the branch was taken to check on the values of 


Z and FIP. 


Listing of Parameters: All paramcters are in COMMON, 

Values taken from and stored in COMMON - the values and 
derivatives are represented by the three arrays STVR, CSTVR, 
and STVRD. The array STVRD corresponds to Z through TEMP 
in COMMON, CSTVR corresponds to CZ through CTEMP, and 
STVRD corresponds to ZDO through TEMDO. T and DT are also 


used, 


1 Program Name: SUBROUTINE HEUN 


2. Program Description: Subroutine HEUN is the numerical 
integration routine, The Runge-Kutta 4th order method is 
used. Very briefly, the method is as follows: we are seeking 
to solve the equation dy/dx = f(x,y). Given an (x Yn we 


= + + + + 
i, (k, + 2k, + 2k k,)/6 


calculate: y 9 3 


ntl 


{ where: k, = hf(x,, ya 


if = + 
ky hf(x +h/2, Ye k,/2) 
ky = hf (x, + h/2,y,, + k,/2) 


= + + 
k, hf (x h; Ve kg) 


f _and_ h = independent axis increment (dx), 


= a Program Logic: Before subroutine CALVR is called each time, 
! a check is made on the values of Z and FIP. If Z is greater than 
[ » 900 Z is set to 1,000 if FIP is negative, FIP is set 'to zero. 
: The local variable IJAX is used to return control to ‘He gaint in 
ne routine where the branch was taken to check on the values of 
~ and FIP. 


4, _. sting of Parameters: All parameters are in COMMON. 
Values taken from and stored in COMMON - the values and 
derivatives are represented by the three arrays STVR, CSTVR, 
¢.d STVRD. The array STVRD corresponds to Z through TEMP 
ij: COMMON, CSTVR corresponds to CZ through CTEMP, and 
S?VRD corresponds to ZDO through TEMDO. T and DT are also 


used. 


sa Peegre ni Tug the 
a. Source program - 84 ~~ statements (including comments). 
b. Object program - 354 machine words. 

6. Programs Called By This Program: Subroutine CALVR. 

at) Programs Calling This Program: EXEC mainline. 

8. Local Variables: 

IJAX - takes on the values once (1) through 4 to 
control the 4 times subroutine CALVR is 
called. 

DD - set to DI/2 for use in solving for k, and 
ky. 

| 9. Piinted Output: None 
eB 
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Program Name: SUBROUTINE CALVR 


Program Description: Subroutine CALVR updates the values 

of all the working variables in the program. This is accom- 

plished by calling the following subroutines: CON1, CON2, QDEV, 
FLUXG, FLUXM, HULL, PROP, QICE (if the ship has encountered 
ice), and STATE. Interm:xed between calls to the above subroutines, 
subroutine CALVR does some calculations and adjusting of its own. 
This subroutine is called by subroutine HEUN (the numerical inte- 


gration routine) each time an updated set of values is desired. 


Program Logic: Subroutine QICE is called only after the ship 

has encountered ice which is controlled by the variable TICE - 

the relative time since the beginning of the run before the ship 
encounters ice. The variable IPLUG (equivalenced to ISPAR(2)) 
is initially set to one (1) so that the initial values of T and FMI will 
be saved for use in the calculation of the updated values of FMI. 


Thereafter IPLUG is set to two(2) to immediately calculate a new 
os ge . “4 é 


a . 
* 


value of FMI using the saved initial values. 


Listing of Parameters: All parameters are in COMMON. Values 
taken from CCMiMON-XNGR, NXPM, NXGRD, XNPMD, QDE, QBRG, 
CWM, XNPM, XGEN, TRNFG, FGI, TRNSF, AJ, CFG, CFM, XNP, RG, 
IPLUG, T, FMI, TAUA, TAUB, FMIRV, THETA, C1, ¥K1, FIMAX, 
FIP, CSHFR, EB, WDED, VS, QPI, TDED, XPROI, XKK, CT, and 
TICE. 

Results stored in COMMON - XNE, XNEDO, QDE, AMPTR, EG, EC, 
QG, EB,QM, FMI, AIHAT, QF RM, VA,QP, TP, and QI. 


-164- 


Program Length: 


a. Source program - 110 statements (ineluding comments). 


b, Ob: tt program - 310 machine words. 


\ 


QDEV, FLUXG, FLUXM, HULL, PROP, QICE, and STATE. 


ts Programs Calling This Program: Subroutines HEUN and CNTRL. 
i] 
8, Local Variables: 
TTN - contains initial value of T. 
1 FMIIC - contains initial value of FMI 


ABLE, BAKER- set to exponential expressions used in 


| the calculation of FMI. 


3. Printed Output: None 


Oey ome 


ioe) 
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Program Name: SUBROUTINE CON1 .- 


Program Description: Subroutine CON1 computes the actual 


bridge lever position (@) and bridge lever time rate (de/dt) which 


depends on time (t) and control values Ty and dpe 


Program Logic: The control values Ty and Tp determine the 
speed of the lever in moving from the initial position (80) to 
the final or ordered position (9*). Thelever moves from 6 
to 6, in i seconds and from @ 


1 1 
is made to determine where the present value of t lies, after- 


‘o oe in qt seconds. A check 
which @ and d0/dt are computed accordingly. 

isting of Parameters: All parameters are located in COMMON, 
‘Values taken from COMMON - T, TAU1, TAU2, THET@, THET1, . 
and THEIC 


Results stored in COMMON - THETA and THETD. 


Program Length: 


a. Source program - 45 statements (including comments). 


b. Object program - 82 machine words. 
Programs Called By This Program: None 

Programs Calling This Program: Subroutine CALVR. 
Local Variables: None 


Printed Output: None 
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i 
| 1, Program Name; SUBROUTINE CON2 | 
| | | 
| 2. Program Description: Subroutine CON2 computes the ordered generator 
[ field voit: - ta ordered diesel speed Se and ordered diesel 
speed rate dn / a0. The values calculated are dependent on the 

f values of the actual bridge lever position (0) and bridge lever time 

| rate (dé/dt) computed in subroutine CON1. 

; | 5. Program Logic: A maximum value of Cre is calculated from the 

Ge generator field resistance Re) and variable resistor in the generator 
field (R,). If the value of 6 is below 63. 0 degrees, the diesel is 
still idling and the value of dn pf at will be zero, If 6 is greater than 
63.0 degrees, a check is made to determine if a full (@=135, 0 degrees) 
or maximum (@=165,. 0 degrees) command was given. The sign of . 

| Cbg is determined to be negative if 6 is negative and positive if 8 
i is positive. 

1 | 
i 4, Listing of Parameters: All parameters are located in COMMON, 

Values taken from COMMON RX, RFG, THETA, THET@ and THETD. | 
if Results stored in COMMON - EFG,XNGR, and XNGRD. 

| [ 5. Program Length: 

| | a. Source program - 56 statements (including comments). 

b. Object program - 178 machine words. 
: j 
; 6, Programs Called By This Program: None 


7. Programs Calling This Program: Subroutine CALVR. 


8, Local Variables: 


EFGM - real variable containing the computed maximum 
value of e, . 
fg 


THET - real value containing the absolute value of 0. 


Printed Output: None. 


1; Program Name: SUBROUTINE FLUXG 
| 
2. Program Description: Subroutine FILLUXG calculates the | 


generator flux (C a” generator field inductance (L, g”’ generator | 
time constant ( r, Q: mutual inductance between ee separate 

and in series faite (XM), and the inductance of the generator in 

series field (XLASF) from the value of ampere-turns (AMPTR), 2 

table of field flux (CFMA), values of control field turns \TRNFG), 

series field tunrs ([RNSF), field resistances (RX, RG) and flux con- 


version constant (ZGP). 


3. Program Logic: The generator field flux is taken trom the input 


array (CFGA) with the value of the generator spacing (SPGEN) being 


the independent axis. Since the ind ependent axis runs from zero to 


12 it is necessary to truncate the independent axis index (KG) and then 


compute the enclosing indices by adding two and one to the value. 


4, Listing of Parameters: All parameters are in COMMON. 
Values taken from COMMON-AMPTR, SPGEN, CFGA, TRNFG 
ZGP, RX, RG, and TRNSF. 


Results stored in COMMON-CFG, XLFG, TAUFG, XM, and XLASF. 


5. Program Length: 


a. Source program - 52 statements (including comments). 


b. Object program - 148 machine words. 


8. Programs Called By This Program: None 


He Programs Calling This Program: Subroutine:CALVR ‘a 


8. Local Variables: 


PG real variable containing the actual value of the 


independent axis index. 
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KG - integer variable containing the truncated value 
of FG. 
DCGDI - real variable containing the interval that the value 


of FG falls in. : ‘ 


9, Printed Output: None 


1, Program Name: SUBROUTINE FLUXM 


2, Program Description: Subroutine FLUXM calculctes the motor 
[ ; field flux (Cen the motor fieid time constant ( Tem)’ and the 
motor field inductance (Len) from the value of field current (FMI), 


[ a table of field flux (CF MA), values of field turns (TRNFM), field 


ee resistance (RIM) and flux conversion constant (ZMP). | 


3. “Program Logic: The motor field flux is taken from the input array 
(CFMA) with the value of the motor spacing (SPMOT) being the | 
independent axis. Since the independent axis runs from zero to 


| 12, it is necessary to truncate the independent axis index (KM) 


the value. 


| 
and then compute the enclosing indices by adding two and one to | 
4, Listing of Parameters: All parameters are in COMMON. 


Values taken from COMMON-FMI,SPMOT, DFMA, TRNFM, ZMP, 
| and RFM. 
Results stored in COMMON-DFM, XLIFM, and TAUFM. 


| 5. Program Length: 


[ a. Source program - 52 statements (including comments). 

“ b. Object program - 112 machine words. 

| | 

; 6. Programs Called By This Program: None 

| Ga Programs Calling This Program: Subroutine CALVR. 

| 8. Local Variables: A 
| FM - ' real variable containing actual value of 


independent axis index. 


KM- integer value containing truncated value of FM 
DCMDI - real value containing the interval th:i th« value of 


IM falls in. 


9, Printed Output: None 


b—aeeg 


otes 


Program Description: Subroutine QDEV calculates the developed 
diesel torque Qi? using the developed torguc constant (K,), the 
servo stroke position (Z), the number of gencrators (XGEN), and 


the angular velocity of the prime mover Na 


Program J.ogic: The maximum developed diesel torque (Qde sax) 


is computed from the angular velocity of the prime mover (Nn 
A diesel engine friction torque is computed using afriction torque - 
ere curve with a Y-intercept indicated by the value of QBRG and a 
slope indicated by CWM. ‘This valuc is added to the old value of 


Que aN ae to arrive at an augmented Qye If the value of Qge is 


max’ 


greater than Qye ase Que is set equal to Qdemax ; 


Listing of Parameters: All parameters are in COMMON, 
Values taken from COMMON- XK6,Z, XNPM, XGEN, QBRG, and CWM. 
Results stored in COMMON - QDI: 


— —— ————. __ 


a. Source program - 53 statements (including comments). 


b. Object program - 138 machine words. 


a a nD NS 


Programs Calling This Program: Subroutine CALVR 


Local Variables: 


QDEMX - rea] variable containing the maximum 


developed diesel] torque. 


Printed Output: None 


Program Name: SUBROUTINE HULL 


Program Description: Subroutine HULL computes the ship's 


resistance (R), wake deduction (WDED), and thrust deduction 


(TDED) from the ship's speed (V.), ice resistance factor (CICE) 
and input tables of effective horsepower (EHP), wake deduction 


(WADED) and thrust deduction (THDED). 


Program Logic: If the control value ICNTR implies ahead 
full in ice (equal to 4), the ice resistance factor (CICE) is taken 
into account in computing the resistance of the hull. If the ship's 
speed is negative, the hull resistance is taken as 1.4 times the 
computed value and the wake and thrust deductions are taken as 


1.0 and 0. 9 respectively. 


Since the independent axis of the three tables (EHP, WADED, and 
THDED) runs from zero to 19.0, it is necessary to add two and 
one to the truncated index (KV) in order to arrive at the interval 


the actual] index (V) falls between. 


Listing of Parameters: All parameters are in COMMON. 

Values taken from COMMON-VS, EHP, WADED, THDED, ICNTR, 
and CICE. | 

Results stored in COMMON - WDED, TDED, and R. 


Program Length: 


a. Source program - 52 statements (including comments). | a 


b. Object program - 152 machine words. 


Programs Called By This Program: None 


ES nn ee 


Programs Calling ‘This Program 


Local Variables: 


is 
KV - 


REHP - 


- Subroutine CAI.VR 


real variable containing the absolute value 


of the ship's speed in knots. 

integer variable coniaining truncated value 
of V. 

real variable containing the effective horsepowcr 


used in calculating NR. 


Printed Output: None 


Program Name: SUBROUTINE PROP 


propeller data tables to compute the ideal propeller torque 
(QP) and thrust per propeller (XPROI) given the propeller 


angular velocity fy advance velocity of the propeller (Va). 


Program logic: The values of the propeller thrust and torque 
coefficients are computed from the tables by a complex but 
efficient means. Depending on the signs of VA and XN, a value 
of the variable KB is set (sec Local Variables). This value is used 
in computing the table subscript and also to set the sign of the 
interpolation fraction (PP) to the negative. A value of XJ or 1/XJ 
is computed, The integer variable KD is set to the truncated 
value of 10 times | xu fe This will determine the relative position 
of 1 in its quadrant, The actual value of Jin the table is computed 
using KD and KB, The interpolation fraction (PP) is computed 
and its sign adjusted by using KB to be negative, The coefficients 


can now be computed from the tables. | 


Listing of Parameters: All parameters are in COMMON. 
Values taken from COMMON -RO,D,XN, VA,B, BB, i, and EE. 
Results stored in COMMON - QPI, and XPROI. | 


Program Length: 


a. Source program - 106 statements (including comments). 


b. Object program - 430 machine words. 


= er pe One ree arene nee ee 


a. Programs Calling This Program: Subroutine CAIL.VR 


8. Local Variables: 


KB - takes on the values zero (0) through 3 


depending on the signs of XN and VA as 


follows. 


I - index to data arrays (B, BB, I, and EE). 

Se ee equal equal to zero (0) if using J, 1 if using © 
1/J 

XJ - equal to VA/(XN- D) if result within range 


/XJ/= 1.5, otherwise equal to XN: D/VA 


KD - equal to 10- /XJ/ and used in computing I 
PE negative fraction used in interp: lation 

XKQ - propeller torque coefficient 

tel = propeller thrust coefficient 

XKQD - modified propeller torque coefficient 

XKTD - modified propeller thrust coefficient 


3 9, Printed Output: None 


| 
2. Program Description: Subroutine QICE computes the ice torque 
(Q)) of the propeller, lf necessary, subroutine RBMAX is called 
to determine a new value of the maximum depth of cut in feet that 
the propeller wil] make in the ice, The value of Q; is taken as the 
lesser of crushing torque (Q.cp) and combined crushing and shearing 
torque Q apg) except that Q, will never be less than the minimum ice 
| torque (Qs vat): 
3, Program Logic: ‘The variable JELO is initially set to one (1) 
in the EXEC mainline so that subroutine RBMAX will be called 
immediatcly the first time this subroutine is executed. There- 
after, subroutine RBMAX is called only when a new value of the 
depth of cut is needed (when the angular position of the propeller is 


g&.eater than the cut angle, ps and greater than or cqual to the 


4. Listing of Parameters: All parameters are located in COMMON, 


Values taken from COMMON-JELO, RBMX, D,1ZB, PI, SIE, XNP, 
PDR, DAR, DCR, FO, EK, DCT, DSH, and VS. 


| 

| ur gle between propeller baldes. 
Results stored in COMMON-JEJ.O and QI, 
t 


5, Program Length: 
| . 
a. Source program - 75 statements (including comments), 
b, Object program - 352 machine words. 


6. Progra. * Called By This Program: Subroutine RBMAX 


Progeams Calling This Program: Subroutine CALVR 


8. 


9, 


Local Variables: 


CSAB2 - 
SNAB2 - 
AB2 - 
its = 


ABB - 


CSBET = 
RB - 
RM - 


ALFAV- 
ALFAI - 
El - 
BAV - 
QICRS - 
QIMIN - 


QICR - 


real variable containing the cosine of a, /2 
real variable containing the sine of a /2 
real variable containing the value of a 
real variable containing the number of propeller 
blades. | 

real variable containing the number of radians 
between prop blades. 

real variable containing the cosine of 6 

real variable containing value of Ry 

although appearing in COMMON, this local 

real variable contains the value of Ra 

real variable containing value of a, 
real variable containing value of Qs 
real variable containing value of ey 
real variable containing value of B ve 
realva le contai:.ing value of QRS 
real variable containing value of Qs MIN 


real variable containing value of Q.cr 


Printed Output: None 


Program Name: SUBROUTINE RBMAX 


oor ae Ske: ey serra. eeroeme + 


muin depth of cut in feet that the propeller blade maker in a block 
of ice, This value is chosen at random from an array of values 
that represent the maximum depth of cut as a percentage of the 


propeller diameter, 


Program Logic; | ¢ ice pattern is made up of two blocks of 


ice ina repeati |: equence, Four variables control the spacing 
and sine of the blocks (BLOC1, BLOC2, DBB12, and DBB21.) 

The value of the variable KELO determines where in the ice 
pattern the ship is at a given instance. The value of KELO is 
initially set to one (1) in the EXEC mainline. When this sub- 
routine is initially executed, the ship had just entered block one and 
KELO is set to two (2). The variable XICE is used to save the 
ship’. distance upon entering the ice to determine how long thc 

ship will remain in a particular block of ice. The variable KELO 
assumes other values depending on whether the ship is between blocks 
one and two, in block two , or between blocks two and one, The 
first 10 values of the RBA array are used when in block one, the 


remaining 15 values when in block two. 


Listing of Parameters: All parameters are in COMMON. 
Values taken from COMMON-KELO , CXS, BLOC1,RBA,D, 
DBB21, BLOC2, and DRB12. 

Results stored in COMMON - CSIE, SIE, RBMX, and KELO. 


Program Length: 


a. Source program - © 87 statements (including’ comments). 
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b, Object program - 178 machine words. 


Proprams Called By This Program: None 


Programs Calling This Program: Subroutine QICE 


J.ocal Variables: 


oa 


real variable used to save the present value 
of the ship's distance cach time a block is 


entered or exited from. 


integer variable used to reference the RBA 
array in determining the maximum ‘Jepth of 
cut; updated by one (1) each time a new value is 


needed. 


Be wee 


ti | 


Program Name: SUBROUTINE STATE 


Program Description: Subroutine STATE calculates the values 

of all the derivatives using the values calculated in the previous 
subroutines. The derivatives calculated are: 1) dz/dt -servo 
stroke velocity; 2) deoal servo stroke acceleration; 3) dn n/t - 
time rate of prime mover speed; 4) di, _/at - time rate of generator 
field current; 4) di',/dt - time rate of motor field adjusting current; 
6) dn /at -propeller acceleration; 7) dv ,/dt - ship's acceleration; 
8) dx _/dt - ship's velocity; 9)SINDO - peopeller angular velocity; 


and 10) di, /dt - time rate of armature circuit current. 


Program Logic: If the value of the actual bridge lever position 
(6) is negative or if the bus voltage (e)) is less than or equal to the 


constant C,, then the time rate of the motor field adjusting current is 


4? 
set to zero. 


Listing of Parameters: All parameters are located in COMMON. 
Values taken from COMMON, Y, ALP1, ALPZ, ALP3, ALPD, XNEDO, 
XNE, Z,QDE, QG, PI,XJDG , EFG, RFG, RX, FGI, EG, EC, RA, Al, 
XLA, XGEN, XLASF, XM, XLFG, THETA, EB, C1, AIHAT, XLFM, QM, 
QFRM, QP, QI, XJMP, TP, R, XMASS, VS and XNP. 

Results stored in COMMON -ZDO , YDO, XNPMD,AIDO, FGIDO, 
FIPDO, XNPDO, VSDO, XS)D0, and SIEDO: 


Program Length: 


a. Source program - 66 statements (including comments). 


b. Object program - 252 machine words. 


Ey 


beac 3 


6, Programs Called By This Program: None 


7. Program Calling This Program: Subreutine CALVR. 


ee ee 
w 


Local Variables: 


VOLT1- real variable used to contain the value of 


fp - (Re + R,) ee 


VOLT2- real variable used to contain the value of 
ear e. - Roi, 
XLD - real variable used to contain the value of 
Lo +1,.5- XLASF/XGEN 


} 
et Fae em kA 
rs 


XLM - real variable used to contain the value of 
2,25 XM /XGEN 


DENOM - real variable used to contain the value of 
1,5 Leg (XLD-XLM) 


$n Si cna aoe aed AE SO 


iA Printed Output: None 


Byte 
chs 


5 
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Program Name: SUBROUTINE OUT 


Program Description: Subroutine OUT prints out the program 


variables in question on the II3M 1132 printer. 


Program Logic: Before printing the variables in question, the 
following values are changed to different units: the value of the 
angular velocity of the prime mover Oa is changed from 
revolutions per second to revolutions per minute as is the 
propeller angular velocity (ns the ship's spced vo) is con- 
verted from feet per second, to knots; the total propeller thrust 
(Tr) and ship resistance (R) are changed from pounds to long 
tons. T'nhe shaft horsepower is also computed from the motor 


torque (Q,,) and a 


Se 


Values taken from COMMON - CXNPM,CXNP, CVS, TP,R, QM, PI, 
XNP, CT, CFGI, QDE, QG, EB, FM1, EC, AI, and QP. 
Results stored in COMMON - None 


Program J.ength: 


a. Source program - 49 Statements (including comments), 


b. Object program - 132 machine words. 
Programs Called By This Program: None 


Programs Calling This Program: EXEC mainline. 


~ 


4 


Local Variables: 


-ipd- 


RPMPM - 


RPMP - 


CVK - 


de a 


i te 


POW - 


Printed Output: 


real variable containing angular velocity of 


prime mover in revolutions per minute. ' 


real variable containing propeller angular 


velocity in revolutions per minute 
real variable containing ship's speed in knots. 


real variable containing total propeller thrust in 


long tons. 


real variable containing ship's resistance in 


long tons. 


Real variable containing ship's horsepower. 


This program prints out the 16 variables in 


question on the IBM 1130 printer. The values occupy one line 


in the following format: 


RPvP CvK/ Qde 


/ if 
bxxx, Xxxbxxx. X xoek. xDXXXX. XDXXX. X XXXXXX, Xx 


& ifm ec m 
xxx. xbxxxx. xbxxx. xbxxxx. cheers. x 


/ POW 


ia Qp / ~pr O\ 
bxxxxxx, XDXXXXXXX. X XXXXxX, XDxXXXX. xX XXXXXX, X 


1, Program Na:ne; SUBROUTINES FINIS 


2. Program Description: Subroutine IINIS interrogates data switches 1 


through 7 to achieve the following results: 


Switch 1 on implies punch new initial condition cards, 2 on implies 
plot data, 3 on implies continue present solution, 4 on prevents 
future typing of data switch setting instcuctions. 5 on implies a 

new value of the time step (dt), 6 on implies a new value of the 
print interval, and 7 on implies a new value of the armature circuit 


inductance (L.,). 


3. Program Logic: If a graph is indicated, the GRAPH mainline 
is loaded and executed by use of a CALI. LINK (GRAPH) statement 
in this routine, The terminating time value (TEND) is set to the 


last time value considered before the graph routine is executed. 


4, Listing of Parameters: All parameters are in COMMON. | 
Values taken from COMMON - THETA, XNPMD, Y, Z, FGI, XNPM, : 
FIP, XNP, VS,AI, AIDO,CT and DT. | 


Results stored in COMMON - TEND, DT, NOu,XLA, and ISPAR(1) 


5. Program Length: 


a. Source program - 92 statemerits (including comments). 


b. Object program - 498 machine words, 


NR Ee RI REE sr _ 


6. Programs Called By This Program: The GRAPH mainline may 


be loaded and executed by use of a CALL LINK (GRAPH) statement. 


Programs Calling This Program: XC mainline 


a a eee 


8. Local Variables: KHY1, KEY2, KEY3, KEY4, KEYS, KEY6, KEYT- 


interger variables containing the switch settings for data switches 


1 through 7 respectively. A value of one (1) indicates the switch 


was on, a value of 2 indicates off, 


If switch 4 was off, the following switch setting 


instructions are typed out: 


TURN ON SW1LCH -- 


| ie TO PUNCH NEW INITIAL CONDITION CARDS 
ie TO PLOT DATA 
3 - TO CONTINUE SOLUTION OF CURREN'T PROBLEM - 
| SWITCH OFF STARTS NEW PROBLEM 
4 - TO SUPPRESS FUTURE TYPING OF THESE 1 
INSTRUCTIONS 
5 TO CHANGE TIME STEP 
6 - TO CHANGE PRINT INTERVAL 
7 - TO CHANGE XLA 
bo TURN ON SELECTOR SWITCHES AND PRESS START 


If new initial condition cards are punched, they arealgo printed 


on the IBM 1132 printer. | 


typed out: 


( Ve 
Sone a 
( its 


If new values of either dt, NON, or La are to be entered, any | 


of the following messages, whichever be appropriate, will be 


Mn. 


TIME STEP 
NON 
XLA 
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sigs, Koh ete cas | Se, 


1, 


2. 


plotter scale factor is reduced by a factor of 10 to plot the remainder 


Program Name: Mainline Routine GRAPH 


Program Description: Mainline program GRAPH controls the 
plotting of the data stored on the disk by the EXEC mainline pro- 
gram, The program first requires that the code numbers of the 
variables to be plotted be entered from the console keyboard. A 
request is then made to determine if a grid is to.be piotted by re- 
quiring a code nuraber to be entcred from the console keyboard. If 

a grid is to be plotted, the program requires three lines of graph 
identification to be entered from the keyboard after which the in- 
formation as well as the graph serial number (computed or read in 
by the GENIS mainline) is printed out along with the grid. The pro- 
gram then loops through the array of code numbers to plot each vari- 


able inturn. If the curve extends beyond 20 seconds, the time axis 


of the curve. Upon the conclusion of the plotting for each curve, . 
the pen is positioned for printing out the curve label. When ali the 
requested variables have been plotted, the program loads and exe- 


cutes the GENIS mainline to begin a new problem. 


Program Logic: When the curves extend beyond 20 seconds, it is 
necessary to compute a pseudo-origin for the remainder of the curve 
since the time axis plotter scale factor is reduced by a factor cf 10 from 
the 20 second mark until the end of the graph. The new origin ‘KDARN', 
is computed by multiplying the last time value encountered less than or 
equal to 19.9 seconds by 10 and subtracting the original value from 

the result. Following time values are added to this pseudo-origin to 
arrive at an augmented time value, 'CTT', for plotting the correspond- 


ing value of the variable in question on the curve. 
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Listing of Parameters: All parameters are located in COMMON. 
Values are taken from COMMON-SPARE(4), TEND,SX,SY,CT and VAR. 
Results stored in COMMON-NOVAR, NVAR, IELO,NIP, and JJ. 


Program Length: 


a. Source program - 246 statements (including comments), 


b, Object program - 738 machinc words. 


Programs Called By This Program: Subroutines GRID, VAPLO, DATA, 
and PRINT. ‘The GENIS mainline is loaded and executed by use of a 
CALL LINK (GENIS) statement following the completion of the entire 
graph, 


Programs Calling This Program: This mainline is loaded and executed 
by use of the CALL LINK (GRAPH) statement in subroutine FINIS (of the 
EXEC mainline). 


Local Variables: 


NPLOT - this integer variable, although in COMMON, is local to this 
program and is used to contain the input from the keyboard 
that determines if a grid is to be plotted; a value greater 
than or equal to one (1) implies plot a grid, a value less 


than one (1) implies no grid is to be plotted. 


TH. character rotation value: real variable containing the number 
of radians from horizontal (+X direction) that a line of char- 


acters is to be plotted. 


We integer variable uwsed to control the up ordown position 
of the pen during the plotting of a variable; initially set 


to one (1) so that the pen can be moved to the first value 


-189- 


ae - 


to be plotted in the up-position for a particular variable, 
thereafter assuincs a value of 2 since the pen will re- 


main down in plotting the remainder of the curve. 


XD, YD-_ real variables used to position the pen for plotting of the 
curve label after the curve has been plotted; set to start the 
first character 0,05 inches to the right and below the last 


value plotted. 


The following local variables come into use only if the curves extend 


beyond 20 secends: 


SXX - time axis plotter scale factor: rea] variable set to one-tenth 
of the initial time axis plotter scale factor, 'SX', for plot- 
ting of the curv?s should they extend beyond 20 seconds 


since after 20 seconds the increments become 10 seconds. 


XD - real variable set to 10 times the actual time value for the 


first time value exceeding or equal to 19, 90 seconds. 


XDARN- real variable set to the value of "XD! minus the first time 
- value exceeding or equal to 19, 90 seconds. This value will 
then be used as a pseudo-origin for plotting the remainder of 
the curve since the plotter scale factor has been reduced by 


a factor of 10, 


CTT real variable equal to the actual time value, 'C’'T', plus the 
value of 'XDARN' for each time step in plotting a curve after 


the 20 second mark, 


oe ee ae oe oe 


oul on the console typewriter: 
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Al--1,EB--2, XNPM--3, FMI--4,XNP--5, XKW--6, 
CVK--7,QPI--8, XPROI--9, CXS--10, Qi--11, TEMP--12 
WHICH VARIABLES TO BE PLOTTED (10 ONLY) - 
(«MMMM MM MOD 
DO YOU WANT A GRID (0-NO, 1- YES) 
() 
If a grid is rcquested, three lines of not more than 50 characters 
each are required to be input after which they will be plotted out 
preceeded by the following information: 


SERIAL NUMBER - - xxxxxxxxxx 
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Program Name: SUBROUTINE GRID 


Program Description: Subroutine GRID is called to plot and label 

the time axis which is represented by the +X direction and 

the ordinate axis which is represented by the +Y and -Y directions. In 
plotting the tirne axis, the first 20 seconds are plotted in increments of 
one second. A test is then made to determine if the axis is to extend 
beyond 20 seconds and if so, the additional length is determined. The 
time axis is then plotted from the 20 second mark to the end in incre- 
ments of 10 seconds, ‘The ordinate axis is plotted from +5,0 to -5,0 
with each unit representing one inch on the graph. The ordinate axis 

is labeled with the scale factors (quantity of each variable corresponding 


to one inch in the Y-direction) for the variables to be plotted. 


Program Logic: In labeling the time axis, the local real variables 'XD' 
and 'WYK' are set to control the positioning of cach label so that they 
will be 0,05 inches below and centered on each tick mark. if the time 
axis extends beyond 20 seconds, the additional length is computed as 
the number of tick marks beyond the one at 20 seconds (one less than 


the number of intervals between tick marks). 


Listing of Parameters: Al) parameters are located in COMMON. 


CS 


Values taken from COMMON-SX, PI], TEND, NOVAR, NVAR, and SY. 
Results stored in COMMON-NIP. - 


Program Length: 


a. Source program - 223 statements (including comments). 


b, Object program - 696 machine words. 


Programs Called By This Program: Subroutine VAPLO. 


Programs Calling This Program: GRAPH mainline. 


Local Variables: 


SXX - 


Ti 


XD - 


WYE - 


TEST - 


TEST1 - 


time axis plotier scale facior: real variable initially set to 
'SX‘ for the plotting of the first 20 seconds of the time axis 
and then to one-tenth of 'SX' for the plotting of the remainder 
of the time axis since after 20 seconds the time interval 
becomes 10 seconds over the same actual distance initially 
used for the one (1) second intervals. 

character rotation value: real variable containing the number 
of radians from horizontal (+X direction) that a line of char- 


acters is to be plotted. 


label centering value: real variable used to center the axes 
labels at each tick mark; equal to 0,05 divided by the plotter 


scale factor for the axis in question. 


real variable used initially to place the X-axis labels 0,05 
inches below each tick mark and then to indicate the relative 


position of each label for the +Y and -Y directions. 


real variable initially used to label the time axis and then 
used to start the Y-axis labels 0,77 inches to the left of the 
Y-axis. Finally used to indicate the starting position (to 
the left of the Y-axis) that each curve actual scale factor 


will be listed on the graph. 


real variable equal to the extent of the time axis beyond 20 


seconds, 


real variable equal to one tenth of the extent of the time 
axis beyond 20 seconds plus a constant of one (1) to plot 


out the remaining tick marks. 
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ITEST - integer variable containing the number of tick marks to 


be plotted beyond the one at the 20 second mark. 


I - integer variable used as a DO-loop index in labeling the 
axes. 
MIKE - integer variable used to store in turn the code number 


for each variable to be plotted for use in a computed GO TO 


statement to plot the actual-scale factor for that variable. 


BY¥1 = real value containing the actual scale factor (quantity of 
variable per inch along the Y-axis) for each variable in 


turn as the scale factor is plotted. 


Printed Output: Any of the following 12 statements concerning the 


Y-axis scale factors may be listed on the graph: 


AI (1 IN, = xxxxx.x AMPS) 

EB (1 IN. - xxx. x VOJ.TS) 
NPM (1 IN, = xxx. x RPM) 
FMI (1 IN. - xxx. x AMPS) 

NP (1 IN. = xxx.x RPM) 

KW (1 IN. = xxxx.x KW) 

VK (1 IN. = xx.x KTS.) 

QP (1 IN. = xxxxxx.x FT-LBS) 
TP (1 IN. 
XS (1 IN. = xxxx. x FEET) 

QI (1 IN. = xxxxxx. x FT-LBS) 
TEMP (1 IN. = xxx. x DEG. CENT. ) 


n 


xxxxxx. x LBS) 


The following time axis label is plotted out: 


TIME (SECS) 


a 


o 


Each time axis tick mark is labeled starting at 1.0 and going through 20. 0 
by one's and then through the final value by 10's. Every whole ordinate 


(Y) axis tick mark is labeled from 5,0 through 0.9 to -5. 0. 
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Program Name: SUBROUTINE DATA 


Program Description: Subroutine DATA retrieves from the disk 
file the next time value and next value of the variable currently 
being plotted. The disk file index, 'IELO', is automatically up- 
dated to the next record after the values are retrieved. The time 
value is placed in location 'CT' and the new value of the variable 
being plotted in 'VAR' , both in COMMON. 


Program Logic: The current value of the disk file index (‘IELO') 
is used to reference the record of the disk file that is next to be 
retrieved. Since there is no way to read just one value from the 
disk, the entire record must be retrieved, after which the value 
of 'JELO' is automatically updated by one (1) to point to the next 
record. The code number of the variable currently being plotted 
is found in the word of the variable list array, 'NVAR', currently 
being considered which is indicated by value of the variable list 
array index 'JJ'. A computed GO TO statement is then used to 


return the new value for the variable in question. 


Listing of Parameters: All parameters are located in COMMON. 
Values taken from COMMON-NVAR, JJ, and IELO. 
Results stored in COMMON -IELO,CT, and VAR. 


Program Length: 


a. Source program - 90 statements (including comments). 


b. Object program - 136 machine words. 


Programs Called By This Program: None 


Programs Calling This Program: GRAPH mainline 


Local Variables: 


NOOK - integer variable used to store the code 


number of the variable in question for 


use in a computed GO TO statement. 


The following real variables, although appearing in COMMON, 


are actually local to this program and are used to store the 


value of each variable when the file record is read from the disk: 


Al - armature circuit current 


EB- bus volts 


RPMPM- diesel speed 

FMI - motor field current 
RPMP- propeller speed 
XKW - power 

CVK - vessel speed 

QPI - propeller torque 
XPROI - propeller thrust 
CXS - distance 

QI - ice torque 

TEMP- motor terperature 


Printed Output : None 


ll. Program Name: SUBROUTINE VAPLO 


2s Program Description: Subroutine VAPLO retrieves tne plotter 


scale factor for the variable that is about to be plotted. Tine 


plotter scale factor is the actuai number of inches along the Y-axis 


that corresponds to one unit of the variable in question and is placec 


in 'SY', 


plotted is found in the word of the variable list array, 'NVAR', 


\ 
: 3. Program Logic: The code number of the variable about to be 
: that is currently being examined which is indicated by the value 


of the variable list array indea, 'NIP'. A computed GO TO state- 


ment is then used to return the appropriate plotter scale factor. 


4. Listing of Parameters: All parameters are located in COMMON, 
Values taken from COMMON-NVAR, NIP, SY1, SY2,SY3,SY4, SY5, 
SY6,SY7,SY8,SY9,SY10, SY11, and SY12. 
.- Results stored in COMMON - SY. 
. 5. Program Length: 


a. Source program - 85 statemerts (including comments). 


b. Object program - 100 machine words. 
r 6. Programs Called By This Program: . None 
{ 7. Programs Calling This Program: GRAPH mainline, subroutine GRID. 


‘H 8, iocal Variables: 


NAN - integer variable used to store the code number 


| 5 of the variable in question for use in :a computed 
GO TO statement. 
Printed Output: None 


Program Name: SUBROUTINE PRINT 


Program Description: Subroutine PRINT is executed upon the 
completion of the curve for a particular variable in order to plot the a , 
appropriate curve label for that variable. The curve label is plotted 


starting at the current location of the pen which was set by the GRAPH 


mainline before calling this routine to be immediately following the 


last value plotted for the variable. 


Program Logic: The code number of the variable in question is 
found in the word of the variable list array, 'NVAR', that is currently 
being examined which is indicated by the value of the variable list 
array index 'JJ'. A computed GO TO statement is then used to plot 


out the appropriate curve label. 


Listing of Parameters: All parameters are located in COMMON. 
Values taken from COMMON -NVAR and JJ. 
Results stored in COMMON ~- none. 


Program Length: 
a. Source program - 97 statements (including comments), 
b. Object program - 212 machine words, 


Programs Called By This Program: None 
Programs Calling This Program: GRAPH mainline 


local Variables: 


MIKE - integer variable used to store the code | 
number of the variable in question for use | ' 


in a computed GO TO statement. 


~199- 


Printed Output: This program may plot out any one of the following 


12 curve labels depending on the variable in question: 


--- ARM. CIR. CUR. 
--- BUS VOLTS 
--- DIESEL SPEED 
--- MOTOR FIELD CUR. 
--- PROP SPEED 
--- POWER 

--- VESSEL SPEED 
--- PROP TORQUE 
--- PROP THRUST 
--- DISTANCE 

--- ICE TORQUE 
--- MOTOR TEMP. 
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